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Measurement of inclusive D* and associated dijet 
cross sections in photoproduction at HERA 

ZEUS Collaboration 

Abstract 

Inclusive photoproduction of D*^ mesons has been measured for photon-proton 
centre-of-mass energies in the range 130 < W < 280 GeV and photon virtuality 
Q 2 < lGeV 2 . The data sample used corresponds to an integrated luminosity of 
37 pb _1 . Total and differential cross sections as functions of the D* transverse 
momentum and pseudorapidity are presented in restricted kinematical regions and 
the data are compared with next-to-leading order (NLO) perturbative QCD calcu- 
lations using the "massive charm" and "massless charm" schemes. The measured 
cross sections are generally above the NLO calculations, in particular in the forward 
(proton) direction. The large data sample also allows the study of dijet production 
associated with charm. A significant resolved as well as a direct photon component 
contribute to the cross section. Leading order QCD Monte Carlo calculations indi- 
cate that the resolved contribution arises from a significant charm component in the 
photon. A massive charm NLO parton level calculation yields lower cross sections 
compared to the measured results in a kinematic region where the resolved photon 
contribution is significant. 
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1 Introduction 



In photoproduction processes at HERA, a quasi-real photon (Q 2 ~ 0) is emitted by the 
incoming electron or positron, and interacts with the proton. Such a photon has a hadronic 
component, which can be assigned a partonic structure. At leading order (LO) in QCD, 
two types of process take part in photoproduction: direct photon processes, where the 
photon couples as a point-like particle to a parton from the proton, and resolved photon 
processes, where one of the partons in the photon scatters on a parton in the proton. 
The light quark structure of the photon has been extensively studied in photon-photon 
collisions at e + e~ storage rings whilst there is little information at present on the 
charm content of the photon. HERA jet studies have shown some sensitivity to the gluon 
content of the photon p], || , which is still poorly known. In this paper we present a study 
of charm photoproduction. Here, the direct process is photon-gluon fusion, 7$ — > cc, 
while charm quarks in the parton distributions of the photon and the proton can lead to 
processes of the type eg — *> eg, known as charm flavour excitation. 

The photoproduction of heavy quarks such as charm can be calculated using perturbative 
QCD (pQCD) with a hard scale given by the heavy quark mass or by the high transverse 
momentum of the produced partons. Two types of QCD NLO calculations are avail- 
able for comparison with measurements of charm photoproduction. The massive charm 
approach @] assumes that gluons and light quarks (u,d,s) are the only active partons 
within the proton and the photon, so that charm is only produced dynamically in the 
hard process. In the massless charm approach ||, [| [7] charm is treated as an additional 
active flavour. The massive approach is expected to be superior for p\ ~ m 2 while the 
massless one is expected to describe the data better for p 2 ± ^> m 2 ||, where p± and m c 
are the transverse momentum and mass of the charm quark. In NLO calculations, di- 
rect and resolved components cannot be unambiguously separated. The massless charm 
calculations take into account charm excitation processes and thus predict, for a given 
factorisation scale, a larger resolved component in comparison with the massive calcula- 
tion. Therefore, it is interesting to compare the predictions of these models to data and 
to investigate the sensitivity of the experimental results to the partonic content of the 
photon and specifically to the charm excitation contribution. 

In the analysis described in this paper, charm was tagged by identifying D* ± (2010) mesons 
in the final state via the charged products of their decay. D* mesons are reconstructed 
through the two decay modesQ: 

D* + - D° 7T + - (K- 7T + ) 7T+ (1) 
D* + -> D° 7T+ -> (K~ 7T+ 7T+ 7T-) 7T + . (2) 

The small mass difference AM = M(D*) - M(D°) = 145.42 ± 0.05 MeV @ yields a low 
momentum pion ("soft pion", its) from the D* decay and prominent signals just above 
the threshold of the M(K7T7Ts) — M(Ktc) and M(Kn h-k-ks) — M^KirKir) distributions, 
where the phase space contribution is highly suppressed ||10|| . 

1 In this analysis D* ± (2010) are referred to as D* and the charge conjugated processes are also included. 
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We present measurements of integrated and differential cross sections for D* mesons pro- 
duced in restricted kinematic regions in p® and r] D * . Here r] D is the pseudorapidity of 
the D*, defined as — ln(tan(6 l /2)), where the polar angle 9 is taken with respect to the 
proton beam direction.^ The data sample is larger by more than an order of magnitude 
compared to our previous study []TTJ , which allows an accurate measurement of the dif- 
ferential cross sections in both p® and r] D * and thus a more stringent test of the NLO 
QCD predictions. 

The improved statistics of the D* sample allows, for the first time, the study of dijet 
photoproduction in association with charm. In such events, the fraction x° BS of the 
photon momentum which participates in the dijet production can be measured H. This 
quantity is sensitive to the relative contributions of resolved and direct processes |L2|. In 



LO QCD direct photon events at the parton level have x® — 1, while resolved photon 
events populate low values of x° BS . 



2 Experimental Conditions 

The data presented in this analysis were collected with the ZEUS detector at HERA during 
the 1996 and 1997 running periods, where a positron beam with energy E e =27.5 GeV 
collided with a proton beam with energy E p =820 GeV. The data sample corresponds to 
an integrated luminosity of 36.9 ± 0.5 pb _1 . A detailed description of the ZEUS detector 
can be found in refs. fl3] , KM . Here we present a brief description of the components 
relevant to the present analysis. 

Charged particles are measured by the Central Tracking Detector (CTD) [|HJ which oper- 
ates in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The CTD 
is a drift chamber consisting of 72 cylindrical layers, arranged in 9 superlayers covering 
the polar angle region 15° < 9 < 164°. The transverse momentum resolution for full 
length tracks is (a p± /p ± ) track = 0.005p^ acfc 0.016 (p*_[ acfc in GeV). The CTD was also 
used to establish an interaction vertex for each event. 

Surrounding the solenoid is the uranium-scintillator sampling calorimeter (CAL) fTB] . The 
CAL is hermetic and consists of 5918 cells each read out by two photomultiplier tubes. 
Under test beam conditions, the CAL has an energy resolution of 0.18/a/E for electrons 
and 0.35 fy/E for hadrons (E in GeV). The effects of uranium noise were minimised by 
discarding cells in the inner (electromagnetic) or outer (hadronic) sections if they had en- 
ergy deposits of less than 60 (110) MeV. For cells without energy deposits in neighbouring 
cells this cut was increased to 80 (140) MeV. 

The luminosity was measured from the rate of the bremsstrahlung process e + p — > e +, yp, 
where the photon is measured by a calorimeter [[T7| located at Z — —107 m in the HERA 
tunnel. 

2 We use the standard ZEUS right-handed coordinate system, in which X = Y = Z = 0is the nominal 
interaction point and the positive Z-axis points in the direction of the proton beam (referred to as the 
forward direction). 



2 



The ZEUS detector uses a three level trigger system [14]. At the first level trigger (FLT) 
the calorimeter cells were combined to define regional and global sums which were required 
to exceed various CAL energy thresholds. In addition, at least one CTD track coming 
from the ep interaction region was required. 

At the second level trigger, beam-gas events were rejected by exploiting the excellent 
timing resolution of the calorimeter and by cutting on the quantity Hi(E — p z )- t > 8 GeV, 
where the sum runs over all calorimeter cells and pz is the Z component of the momentum 
vector assigned to each cell of energy E. In addition, events were rejected if the vertex 
determined by the CTD was not compatible with the nominal ep interaction point. 

At the third level trigger (TLT) the full event information was available. Calorimeter 
timing cuts were tightened in order to reject the remaining beam-gas events. At least 
one combination of tracks detected in the CTD was required to be within wide mass 
windows around the nominal values in AM and in M(Ktt) (M(Kinnr)) for reaction (1) 
(reaction (2)). In addition, cuts were made on the transverse momenta of tracks associated 
with these D* candidates and p®* was required to be above 1.8 GeV for reaction (1) and 
above 3.3 GeV for reaction (2). For the measurement of D* in association with jets, an 
alternative trigger strategy is possible at the TLT, based upon the jets themselves. The 
jet reconstruction algorithm used the CAL cell energies and positions to identify jets. 
Events were required to have at least two jets, each of which has a transverse energy 
Excai > ^ GeV and pseudorapidity rf et < 2.5. This strategy is used as a cross-check for 
the results on dijets in association with charm. 



3 Analysis 

3.1 Offline Data Selection 

The event sample was processed using the standard offline ZEUS detector calibration and 
event reconstruction code. To define an inclusive photoproduction sample, the following 
requirements were imposed: 

• A reconstructed vertex with at least three associated tracks. 

• No scattered positron found in the CAL by the algorithm described in ref. ||12| . 
This requirement removes neutral current deep inelastic scattering (DIS) events, 
thereby restricting Q 2 to below ~ 1 GeV 2 . The corresponding median Q 2 in our 
photoproduction sample is estimated from Monte Carlo (MC) simulations to be 
~ 3-10- 4 GeV 2 . 



115 < Wj B < 250 GeV, where W JB = ^Ay JB E p E e . Here W JB and 
Djb = ^i(E — p z )i/2E e are the Jacquet-Blondel [18| estimators of W and y, re- 



spectively, and y is the fraction of the positron beam energy taken by the photon. 
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The value of Wjb was determined from the energy deposits in the uranium calorime- 
ter. The lower Wjb cut rejects events from a region where the acceptance is small 
because of the trigger requirements. The upper cut rejects possible background from 
DIS events in which the scattered positron had not been recognised. A systematic 
shift in the reconstructed values of Wjb with respect to the true W of the event, 
due to energy losses in inactive material in front of the calorimeter and particles lost 
in the beam pipe, was corrected [12, pij, using the MC simulation of the detector 
described in section 4. The centre of mass energy range covered by the photopro- 
duction sample is then 130 < W < 280 GeV, corresponding to 0.19 < y < 0.87. 



3.2 Reconstruction of D* Candidates 

A D* reconstruction algorithm was applied to all selected events. It uses the mass dif- 
ference technique to suppress the high background due to random combinations from 
non-charm events, which have a much higher cross section. Only tracks associated with 
the event vertex and having p l ™ ck > 0.15 GeV and \r] track \ < 1.75 were included in the 
combinations. 

Reconstructed tracks in each event were combined to form D° candidates assuming the 
decay channels ([D or (|2|). For both cases, D° candidates were formed by calculating the 
invariant mass M(Ktt) or M (KTrnn) for combinations having a total charge of zero. No 
particle identification was used, so kaon and pion masses were assigned in turn to each 
particle in the combination. Transverse momenta of p t £ ack > 0.5 GeV were required for 
all tracks of channel ([!]) and for the track taken to be the kaon for channel (Q). Pion 
candidates in the latter channel were required to have p t J ack > 0.3 GeV. An additional 
track, assumed to be the soft pion, 7i s , with a charge opposite to that of the particle 
taken as a kaon, was then added to the D° candidate. The mass difference AM = 
M(KirTts) — M(Kir) for channel (p]) or M(K7nnnrs) ~ M(Knirn) for channel (§) was 
evaluated. The reconstructed D* candidates were required to be in the pseudorapidity 
range — 1.5 < i] D * < 1.5 , for which the CTD acceptance is high. 

To comply with the pj* cut applied at the TLT, we required p^- 71 " 7 ^ > 2 GeV for chan- 
nel ([!]). The number of decay particles in channel (0) is larger: to improve the signal 
to background ratio, we required p ± 7r7r7r7rs > 4 GeV for this channel. Since more combi- 
natorial background exists in the forward direction as well as in the region of low p®* , 
an additional cut, p®* / E^ 10 " > 0.1, was applied to both channels. Here E^ l0 ° is the 
transverse energy outside a cone of 6 = 10° defined with respect to the proton direction. 
This cut, as verified by MC studies, removed a significant fraction of the background 
whilst preserving 99% of the D* signal. 

The AM distributions of channel ([!]) and channel @ for combinations with M(Ktt) or 
M(i^7T7T7r) between 1.80 and 1.92 GeV are shown in Fig.l. Clear peaks at the nominal 
value of M(D*) — M(D°) are evident. MC studies have shown that the contribution of 
other D° decay modes to the AM peak is small and can be neglected. 
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The AM signals were fitted, using a maximum likelihood method, to a sum of a Gaus- 
sian (describing the signal) and a functional form (describing the background shape) of 
A ■ (AM - m n ) B . The mass values obtained were AM = 145.45 ± 0.02 (stat.) MeV for 
channel (|1|) and 145.42 ± 0.05 (stat.)MeV for channel (0), in agreement with the PDG 
value @. The width of the signals were a = 0.68 ± 0.02 MeV and a = 0.72 ± 0.05 MeV, 
respectively, in agreement with our MC simulation. 

To determine the background under the peak for channel (p, combinations in the same 
M(Kir) range, in which both tracks forming the D° candidates have the same charge, with 
its having the opposite charge, were used. These are referred to as wrong charge combi- 
nations. The AM distribution from such combinations is shown as the dashed histogram 
in Fig. la. The inset to Fig. la shows the M(Ktt) distribution from combinations having 
a mass difference in the range 143 < AM < 148 MeV. A D° peak is clearly observed. 
The dashed histogram shows the wrong charge combinations defined above. The excess 
of events with respect to the wrong charge distribution below the D° region originates 
mostly from D° decays involving neutral pions Oj. The number of reconstructed D* 



mesons in channel ([!]) was determined by subtracting the wrong charge distribution after 
normalising it to the right charge distribution in the range 150 < AM < 170 MeV. After 
subtracting the background from the AM distribution of Fig. la, a signal of 3702 ± 136 
D* events was obtained for p®* > 2 GeV. 

Side band subtraction, close to the signal region, was used to estimate the 
background under the AM signal of channel @. The side bands, 1.70 < M^Kinrir) < 
1.80 GeV and 1.92 < M(Ktt'K'k) < 2.02 GeV, were normalised to the region 148 < AM < 
165 MeV (dashed histogram in Fig. lb). This subtraction removed the combinatorial 
background coming from events or tracks in which no D* decaying through this channel 
is produced, and part of the background due to the mass misassignment of the kaon and 
pion candidates with the same charge from the D° decay. The remaining background, 
coming from the mass misassignment, is reproduced in the MC acceptance calculations. 
The inset to Fig. lb shows the M(K7nm) distribution from combinations having a mass 
difference in the range 143 < AM < 148 MeV. A D° peak is clearly observed. The total 
number of D* mesons in channel (2) extracted for p® > 4 GeV from the AM distribution 
with the side band subtraction was 1397 ± 108. 



3.3 Jet Reconstruction 

For the measurement of charmed dijet events, the KTCLUS cluster algorithm |2T| has 



been implemented in its "inclusive" mode [21|. In this algorithm, jets are unambiguously 



defined at the hadron, parton and CAL levels. Using the pt recombination scheme pi] , 
the parameters of the jets are calculated as: E? 1 = X^-E^; rf et = (1/ E^ t ){J^ i E T .r] i )] 
fa> et = (1/ E^ t )(J2i E^fa)- The sums run over all calorimeter cells, hadrons or partons 
belonging to the corresponding jet. Here Et { , t]i and fa are the transverse energy, pseu- 
dorapidity and azimuthal angle. 

For the analysis of charm with associated dijets, events containing a D* meson in chan- 
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nel (1) with > 3 GeV were used. The events were also required to have at least two 

jets with \rji et \ < 2.4 and a reconstructed Ej^ cal > 5 GeV. With this selection, 587 ± 41 
events were found after subtraction of the wrong charge background. In addition, an 
analysis with E^lai > 4 GeV was performed, yielding 971 ± 52 events. The distribution of 
the distance between a D* candidate and the jet closest to it in the rf et - <^ et space shows 
that the measured D* belongs to one of the two jets. In more than 80% of the cases, this 
distance was less than 0.2, which is consistent with the observed hard fragmentation of 
heavy quarks |22 . 



4 Monte Carlo Simulation 



The MC programs PYTHIA 6.1 JEJ and HERWIG 5.9 |2§ were used to model the 
hadronic final states in charm production and to study the efficiency of the cuts used in 
the data selection. Both programs are general purpose generators including a wide range 
of photoproduction processes. 

Large samples of charm events were generated for channels (|l|) and (|2|) using both MC 
programs. Direct and resolved photon events, including charm excitation, were gener- 
ated using as a reference sample the MRSG [^J parametrisation for the proton and 
GRV- G HO |2(J for the photon. These samples have at least ten times the statistics of 
the data, so their contribution to the statistical error is negligible. To check the sensi- 
tivity of the results to the choice of the structure function, the reference samples were 
reweighted to simulate other parton distributions of both the proton and the photon. The 
MC studies showed that, in the kinematic range used here, the results are insensitive to 
contributions from charm excitation in the proton. 

In order to include photoproduced D* mesons originating from b quark events, a sample 
of such events was generated with a ratio to the charm sample proportional to the cross 
section ratio of the two processes used in the MC (~ 1 : 100). Within the kinematic 
range of the inclusive D* analysis, the contribution of b quark production to the D* cross 
section is estimated to be ~ 5%. For the kinematic range of dijets in association with 
charm the corresponding estimate is ~ 10%. 

Events containing at least one D* decaying into channel (|1|) or (|2]) were processed through 
the standard ZEUS detector and trigger simulation programs and through the same event 
reconstruction package used for offline data processing. Tracks were reconstructed both in 
the TLT and the offline simulations. The MC efficiency of the tracking trigger was checked 
using the jet trigger described in section 2 and found to be consistent with the data. 
Satisfactory agreement was observed between the CTD transverse momentum resolution 
in the MC samples and the data. 

An additional sample of events was generated using multiparton interactions (MI) in 
HERWIG [^7| as an attempt to simulate the energy from additional softer scatters ( "un- 
derlying event"). 



6 



5 Measurement of Inclusive D* Cross Sections 



The improved trigger and detector conditions compared to that used for our previous 
results |]T]J allow measurements of the inclusive ep — ► D*X cross sections in a wider 
kinematic region: p® > 2 GeV and —1.5 < i] D < 1.5. The integrated D* cross section 
in the above region for Q 2 < lGeV 2 , 130 < W < 280 GeV was calculated using the 
formula a ep ^o*x = N®* rr /£B, where N®* rr is the acceptance-corrected number of D* , B 
is the combined D* and D° decay branching ratios (0.0262 ± 0.0010 for channel (1) and 
0.051 ± 0.003 for channel (2)) and £= 36.9 ± 0.5 pb _1 is the integrated luminosity. 

In order to obtain N®* rr , a correction factor u>i, defined as the number of generated divided 
by the number of reconstructed D* mesons, was calculated for channel (1) from the MC 
simulation using a three-dimensional grid in the quantities p®* , r] D * and Wjb- The index i 
corresponds to a given grid bin. All D* data candidates in a grid bin were corrected by the 
appropriate uii, yielding N^ rr = Sjtu^A 7 "^*)^ Here {N®* c )i is the number of reconstructed 
D* candidates in bin i. For channel (2) a one dimensional bin-by-bin unfolding procedure 
was used. 

The reference MC used to calculate the acceptance for channel (1) was HERWIG. For 
channel (2) PYTHIA was used, since HERWIG does not reproduce the decay widths 
of resonances which contribute to the Knim final state 0. Results obtained from the 
alternative MC were used in each channel to estimate the systematic uncertainties. 

Table 1 summarises the results for N^* c after background subtraction and the integrated 
cross sections for both decay channels with various p®* cuts. The first error is statistical 
and the second is the combined systematic uncertainty. The overall scale uncertainties 
(±1.4% from the luminosity measurement, and ±3.7% or ±5.7% from the branching 
ratios || of channels (1) or (2) respectively) were not included in the combined systematic 
errors. 

The differential cross sections da /dp®* and da/dr] D were measured using the same pro- 
cedure. The combinatorial background was subtracted bin-by-bin from each distribution 
using the methods described above. The da /dp®* distribution is shown in Fig. 2 for 
— 1.5 < 7] D * < 1.5 for channels ([!]) all and ([|) and listed in Table 2 for channel (1). The 
da/dr] D * distributions for p^ > 2 and 3 GeV for channel ([]]) are shown in Fig. 3(a,b) and 
for p^* > 4 and 6 GeV for both channels in Fig. 3(c,d). In Table 3 the da/drj D values are 
listed for channel (1). 

The results from the two D° decay modes are in good agreement and are consistent with 
our published measurements based on data taken in 1994 [p| . 

5.1 Systematic Uncertainties 

A detailed study of possible sources of systematic uncertainties was carried out for all the 
measured cross sections. The numbers quoted below are for the integrated cross section 
with p 1 ^ > 2 GeV of channel (JT|) , unless stated otherwise. 
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> 2 GeV 


vT > 3 GeV 


pf > 4 GeV 


pf > 6 GeV 




3702 ± 136 


2619 ±82 


1505 ± 50 


410 ±24 


N°l{K<K<KTnx s ) 






1397 ± 108 


411 ±40 



0- dat a(#7r7r s ) [nb] 18.9 ± 1.2 9.17 ± 0.35 +° f Q 4.24 ± 0.16 1^ 0.948 ± 0.061 jftgg 

a d ^(K m <K s ) [nb] 4.22 ± 0.33 ±° f 5 0.991 ± 0.098 t° f e 9 3 

^massive [nb] 13.1 5.43 2.46 0.665 

owicss i [nb] 25.3 8.50 3.37 0.739 

tWiess j [nb] 17.4 5^83 2M O520 

Table 1: Number of reconstructed D* mesons after background subtraction and integrated 
cross sections, a ep ^ D *x, for Q 2 < 1 GeV 2 , 130 < W < 280 GeV, -1.5 < rf* < 1.5 and 
various p®* cuts. Predictions of the NLO QCD calculations are given for the reference pa- 
rameters and parton density functions (see section 6). The first error is statistical and the 
second is systematic. Overall normalisation uncertainties due to luminosity measurement 
(±1.4%) and to D* and D° decay branching ratios (±3.7% for channel (1) and ±5.7% for 
channel (2)) are not included in the systematic errors. 



Uncertainties originating from the modelling of the MC simulation were estimated 
from the difference in the cross sections obtained with the two event generators 
PYTHIA and HERWIG. They are negligible for the cross section with pY" > 2 GeV; 
however they vary between —2.2% and —4.5% for the higher p^ cuts. 

To estimate the uncertainties in the tracking procedure, the track selection cuts 
were varied by ±10% from the nominal values (section 3.2). The resulting combined 
uncertainty in the cross section is Changing the p®* / ' E 6 ^ 100 cut by the same 

amount yields an uncertainty of 1^4% ■ 

The MC simulation was found to reproduce the absolute energy scale of the CAL to 



within ±3% f28|. A shift of ±3% due to the CAL energy scale uncertainty produces 
a variation of in the cross section. The dominant source of this uncertainty 

is due to the acceptance of the CAL energy thresholds in the FLT (section 2). 
An additional uncertainty due to a small mismatch between data and MC in the 
observed CAL energy distribution amounts to 



-1.5a/ 
-1.2 A)- 

Uncertainties in the background estimation of 1q!4% were obtained by varying the 
AM and M(D°) mass windows and the normalisation region (section 3.2). 

The uncertainty from correcting Wjb to the true W, determined by moving the 
Wjb boundary values by the estimated resolution of ±7%, was negligible for the 
cross section with p^* > 2 GeV. For higher p^* cuts the uncertainty varies between 
-1.7% and ±1.5%. 



Reweighting the reference MC samples to other parton density parametrisations |29 
for the proton (MRSA', GRV94HO, CTEQ3M) gave a variation of ±?;g% in the cross 
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pi* (range) GeV da/dpf (nb/GeV) 



9 4^8 


( 
\ 


9 


_ Q 

o 


) 




+1.55 
-0.56 


3.464 


( 


3 


- 4 


) 


4.94 ±0.31 


+0.30 
-0.31 


4.469 


( 


4 


- 5 


) 


2.22 ±0.13 


+0.10 
-0.11 


5.470 


( 


5 


- 6 


) 


1.076 ±0.073 


+0.071 
-0.043 


6.902 


( 


6 


- 8 


) 


0.328 ±0.024 


+0.020 
-0.013 


9.672 


( 


8 


-12 


) 


0.067 ±0.008 


+0.004 
-0.006 



Table 2: The differential cross section da /dp® for channel (1) as function of p±* for 
the kinematic region of Fig.2. The p®* points are given at the positions of the average 
values of an exponential fit in each bin. The p±* range is given in brackets. The first 
error is statistical and the second is systematic. Overall normalisation uncertainties due 
to luminosity measurement (±1.4%) and to D* and D° decay branching ratios (±3.7%) 
are not included in the systematic errors. 



rf* range 


da/dr] D * 


(nb) 




pf > 2 GeV 


pf > 3 GeV 



-1.5, 


-1.0) 


8.89 ±0.81 


+0.89 
-0.33 


2.96 ±0.23 


+0.24 
-0.26 


-1.0, 


-0.5) 


8.16 ±0.80 


+0.70 
-0.48 


4.17 ±0.29 


+0.15 
-0.26 


-0.5, 


0.0) 


7.61 ±0.79 


+0.97 
-0.48 


3.88 ±0.28 


+0.25 
-0.18 


0.0, 


0.5) 


5.23 ±0.99 


+1.09 
-0.55 


2.93 ±0.28 


+0.18 
-0.25 


0.5, 


1.0) 


3.21 ± 1.02 


+0.86 
-0.66 


2.11 ±0.28 


+0.19 
-0.21 


1.0, 


1.5) 


4.65 ± 1.40 


+1.49 
-1.06 


2.32 ±0.34 


+0.46 
-0.37 



pf > 4 GeV 



6 GeV 



-1.5,-1.0 
-1.0,-0.5 
-0.5, 0.0 
0.0, 0.5 
0.5, 1.0 
1.0, 1.5 



1.021 ± 0.096 ±S."oo5 
1.641 ±0.137 IS 
1.877 ±0.143 ±S:ii2 
1.662 ± 0.137 tomr 
1.090 ±0.1 17 IS 
1.186 ±0.142 ±S:i3T 



052 ±0 
331 ±0 
460 ±0 
398 ±0 
374 ±0 
242 ±0 



n c- 9 +0.048 
uoz -0.037 
n47 +0.054 
-0.054 

uui -0.039 
n cr fi +0.032 
uo ° -0.028 
n4 r +0.029 
-0.036 
nfi9 +0.030 
uoz -0.043 



Table 3: The differential cross sections da/di] D * for channel (1) as function of r] D * for 
the kinematic regions of Fig. 3. The i] D * range is given in brackets. The quoted cross 
sections correspond to the centres of the corresponding bins. The first error is statistical 
and the second is systematic. Overall normalisation uncertainties due to luminosity mea- 
surement (±1.4%)and to D* and D° decay branching ratios (±3.7%) are not included in 
the systematic errors. 
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section. Since the photon structure is not well known, we used several parton density 
parametrisations (LAC-G1, ACFGP, GS-G HO) and in addition we allowed a ±10% 
variation of the ratio of resolved to direct photon contributions with respect to the 
reference structure function. The largest resulting uncertainty in the cross section 
was la6%- 

All contributions to the systematic uncertainties, except the overall scale uncertainties, 
were added in quadrature. The combined systematic uncertainties in the cross sections 
are given in Table 1. For the differential cross sections the systematic errors were added 
in quadrature to the statistical and are indicated in Figs. 2-4 by the outer error bars. In 
Tables 2 and 3 both types of errors are given separately. 



6 Comparison with NLO QCD Calculations 



6.1 Massive Charm Scheme 



Full NLO calculations in the massive charm scheme of total and differential cross sections 
for heavy quark production in the HERA kinematic region have been published in ref . [|J . 
The computation was done as in ref. 



Tl|] for 7£> — > ccX |30 



ep — > ccX cross section with the appropriate flux factors [19 
fragmenting into a D* + as measured by the OPAL collaboration [[31 
was used to produce total and differential D 
regions of our measurements. 



and then converted to 
The fraction of c quarks 
0.222±0.014±0.014, 



cross sections in the restricted kinematic 



The calculation used the MRSG [03] and GRV-G HO 0] parton density parametri- 
sations for the proton and photon, respectively. The renormalisation scale used was 
fi R = m±_ = \Jm1+ p\ (m c = 1.5 GeV) and the factorisation scales of the photon and 
proton structure functions were set to fiF = 2m_i_. The charm fragmentation into D* was 

performed using the Peterson function f(z) oc z (1 — 1/z — e/(l — z)Y 



. Here z is 

the fraction of the charm quark momentum taken by the D* and e is a free parameter. 

The NLO cross sections obtained for the same kinematic regions as the data are listed 
in Table 1 for e=0.02 and shown in Figs. 2 and 3 for e=0.02 (dashed lines) and e=0.06 



(dash-dotted lines). The value e=0.06 is based on ref. [3!| and used by Frixione et al. 
in ref. 0, while e=0.02 is suggested by recent fits to e + e" data 0. The predicted cross 
sections are considerably lower than those measured. The dotted line, which corresponds 
to the extreme choice /i# = 0.5m± and m c = 1.2 GeV, is still below the data at the 
high rj D regions. The calculated shapes in both the p^ and r) D * distributions are also 
inconsistent with the data. 

The result of applying an effective intrinsic transverse momentum, kx-, to the incoming 



partons in the massive charm scheme |30| is relatively small. The predicted cross sections 
increase by about 10% with (k T ) = 1 GeV 2 , mostly at low px and in the backward 
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direction. In a semi-hard approach |33| this effect was calculated according to the BFKL 



evolution ||34||. Recently LO predictions using this approach have become available [35 



The predicted cross sections for our kinematic range are close to the data in absolute 
value but do not match the shape of the rf* distribution. 



6.2 Massless Charm Scheme 

A second type of NLO calculation || [| [7[, the massless charm scheme, assumes charm to 
be an active flavour in both the proton and the photon. The two massless charm calcu- 
lations factorise the perturbative and non-perturbative components of the fragmentation 
differently and fit the latter part to the Peterson function ||22|| , using recent e + e~ data on 



D* production to extract the e parameter. The fitted values obtained by the two calcula- 
tions in their specific factorisation schemes are e=0.116 0] and e=0.02 0]. Similar cross 
sections are obtained in each of the massless charm calculations by fitting fragmentation 
functions other than the Peterson one to the e + e~ data. These predictions are expected 
not to be reliable when the minimum p^ cut is as low as 2 GeV. 

The cross sections predicted with these calculations [||, [7| for the kinematic region of our 
measurement are listed in Table 1 and shown as full lines in Figs. 2 and 3. The parton 



density parametrisations used were CTEQ4M for the proton and GRV-G HO [£6[ for 
the photon. The renormalisation and factorisation scales as well as the values of m c are 
the same as in the calculation of the massive charm approach. 

The predictions of the two massless charm models give similar shapes of the differential 
cross sections (Figs. 2 and 3), but disagree with each other in absolute magnitude by 
~ 40%. The cross sections obtained by these predictions are mostly below the data. 
In particular the data are above the NLO expectations in the forward direction. The 
contribution of D* produced from bb in our kinematic region, not included in the NLO 
curves, is predicted [|] to be below 5%, in agreement with our MC estimation (section 4). 
This fraction is found from the MC studies to be slightly higher in the forward region, 
where it is up to 7%. 



Using the MRSG f25| parton density parametrisation of the proton has no significant 
effect on the predictions. In contrast, the calculations depend on the parton density 
parametrisations of the photon and in particular its charm content. In order to check the 
sensitivity of the da/drj D * data to the parton density parametrisation of the photon, we 
compare the results for p±* > 3 GeV and p® > 4 GeV in Fig. 4 with the two NLO massless 
charm predictions [7] obtained with the photon parton density parametrisations GRV- 
G HO [p6fl , GS-G HO |37j and AFG [j38|l . The differences between the various photon 



parton densities are at the 20% level or less in the integrated cross sections, but in the 
differential cross sections considerable differences in shape are observed. For the massless 
charm scheme of ref. |J, the GS-G HO curves |37]] are closest to the data. However, 
in the GS-G HO parton density function used for this calculation, charm and u-quarks 
contribute equally. 
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7 Measurement of D* Dijet Cross Sections 

Given the discrepancies observed between data and NLO predictions in the inclusive D* 
measurements, it is of interest to study the kinematics of charm production in more detail. 
The measurement of jets in the final state allows the kinematics of the hard scattering 
process to be reconstructed. In order to compare the measurement with QCD calculations 



at any order, we define 39 



where rf et is the jet pseudorapidity, y is estimated by i/jb, and the jets in the sum are 
the two highest E T et jets within the accepted rf et range. The variable x° BS is the fraction 
of the photon momentum contributing to the production of the two jets with the highest 
E T et . In measurements, as well as in MC simulations and higher order calculations, direct 
and resolved samples can be separated by a cut on x^ BS . In this analysis we define a 
direct (resolved) photon process by the selection x® BS > 0.75(< 0.75). 

Fig. 5 shows the uncorrected transverse energy flow, (l/Nje^dEr/dAr), around the jet 
axis ("jet profile") as a function of Ar] = r] CELL — rj 3et , the distance in r] of the CAL cell 
from the jet axis for the sample of dijet events associated with a D* (section 3.3) with 
^Tcai > ^ GeV. As for the inclusive D* analysis (section 3.2), wrong charge combinations 
were used to subtract background from the Et flow in the D* signal region. In order 
to reduce the uncertainties due to the background subtraction procedure, a narrower 
D* region was used in the jet profile plots: 1.82 < M(Kn) < 1.90 GeV and 0.144 < 
AM < 0.147 GeV. The jet sample is divided into three regions of r] jet : —2.4 < r] jet < 0.0, 
0.0 < f] 3et < 1.0 and 1.0 < f] 3et < 2.4. The distributions are plotted separately for direct 
(x° BS > 0.75) and resolved (x° BS < 0.75) events. The jet profiles are compared to the 
results of the HERWIG MC which includes LO-direct and LO-resolved photon processes 0, 
shown as the full histogram. In inclusive dijet events |1, the MC simulation gives too 
little transverse energy in the forward (positive A^) region for Iow-Et jets, even when 
that simulation includes MI. In contrast, our charm dijet Et flow distributions are in 
reasonable accord with the MC without MI, including the forward region. 

Also shown in Fig. 5 are the jet profiles obtained if only HERWIG LO-direct photon 
events are used (dotted histogram). These profiles have reduced E T flow in the backward 
(negative Arj) region and do not describe the data with x® BS < 0.75, in particular for 
the ranges < r] jet < 1 and 1 < nf et < 2.4. The Et flow in the backward direction is 
consistent with the presence of a remnant from the resolved photon. 

To calculate the cross section da/dx® BS for dijets with an associated D* meson, MC event 
samples have been used to correct the charm dijet data for the efficiencies of the trigger 
and selection cuts and for migrations caused by detector effects. The resolution of the 
kinematic variables was studied by comparing the MC simulated jets reconstructed from 



3 We distinguish between LO-direct and LO-resolved photon contributions using the LO diagrams as 
implemented in the MC simulation. 
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final state particles (hadron jets) with jets reconstructed from the energies measured in 
the calorimeter (detector jets), and by comparing the corrected jjjb with the true y. The 
resolutions obtained are: in E^ 1 — 15%, in rf et ~ 0.1 and in i^ BS ~ 0.06. The correction 
factors are calculated as the ratio N true /N Tec in each x® BS bin, where iVt rue is the number 
of events generated in a bin and iV rcc is the number of events reconstructed in that bin 
after detector simulation and all experimental cuts. 

Differential cross sections in da/dx° BS in the range 130 < W < 280 GeV , Q 2 < 1 GeV 2 
are given for jets with \rf et \ < 2.4, Ej? > 7 GeV, E^ t2 > 6 GeV and at least one D* 
in the range p®* > 3 GeV, — 1.5 < r] D * < 1.5. The asymmetric cut on the hadron level 
E|f* values has been applied in order to avoid a problem associated with a singularity in 
the NLO calculations due to the soft gluons that accompany the jet The increased 
minimum p®* of 3 GeV compared to the inclusive D* analysis (section 5) is due to the fact 
that there is almost no D* signal in the region below this value due to the requirement 
of the dijet cuts. Background subtraction was performed as described in section 3.2 for 
channel (1). 

The da/dx° BS results are shown in Fig. 6 and listed in Table 4. All uncertainties except 
that due to the energy scale have been added in quadrature. The systematic uncertainty 
due to the energy scale is shown in Fig. 6 as the shaded band. The cross section inte- 
grated over x° BS is 1.65 ±0.12 (stat.) toloe ( s V s t-) -o.w (energy scale) nb. Results are also 
presented in Table 4 for the region E J T etl >6 GeV, E 3 T et2 >5 GeV, where the cross section 
integrated over x° BS is 2.57 ± 0.14 (stat.)tom (syst.) toll (energy scale)nb. 

7.1 Systematic Uncertainties 

Sources of systematic uncertainties in the cross section measurements were investigated in 
a similar manner to section 5.1. Additional contributions specific to the D* and associated 
dijet sample for the integrated cross sections in x)r with E 3 T ea >7 GeV, E J T eU >6 GeV 
are: 

• The possible shift in the CAL energy scale was increased to ±5% due to the addi- 
tional uncertainty in the E^l al measurement 0. The variation in the cross section 

is i^ 2 8 2 %. 

• The uncertainty due to shifting the minimum Ej: l cal cut by ±1 GeV, which corre- 
sponds to the jet resolution in this low energy region is estimated to be to!i%- 

• Varying the rf et cut values by ±0.1 yields an uncertainty of t!\\%- 

• Using the HERWIG MC with MI for the acceptance calculations contributes an 
uncertainty of ±1.1%. 

All contributions to the systematic errors, excluding luminosity, branching ratios and 
energy scale uncertainties, were added in quadrature. The final systematic uncertainty in 
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x° BS range 


da/dx° BS {nb) 


da/dx° BS {nb) 




E 3 T etl > 7GeV 
E 3 T et2 > 6 GeV 


E 3 T etl > 6 GeV 
E J T et2 > 5 GeV 



(0.000- 
(0.125- 
(0.250- 
(0.375- 
(0.500- 
(0.625- 
(0.750- 
(0.875- 



0.125 
0.250 
0.375 
0.500 
0.625 
0.750 
0.875 
1.000 



0.32 ± 
1.06 ± 
1.20 ± 
0.98 ± 
1.24 ± 
1.80 ± 
3.70 ± 
2.87 ± 



0.19 
0.30 
0.28 
0.31 
0.27 
0.36 
0.53 
0.37 



+0.14 
-0.26 
+0.17 
-0.22 
+0.17 
-0.36 
+0.26 
-0.26 
+0.33 
-0.24 
+0.48 
-0.20 
+0.61 
-0.65 
+0.36 
-0.33 



+0.00 
-0.06 
+0.10 
-0.13 
+0.17 
-0.14 

+0.15 
-0.11 
+0.18 
-0.12 
+0.24 
-0.19 
+0.54 
-0.37 
+0.23 
-0.18 



0.42 ± 
1.80 ± 

1.64 ± 
1.58 ± 
1.92 ± 
2.97 ± 
6.34 ± 
3.86 ± 



0.21 
0.35 
0.33 
0.38 
0.34 
0.44 
0.65 
0.42 



+0.33 +0.06 
-0.23 -0.08 
+0.53 +0.30 
-0.85 -0.20 
+0.52 +0.24 
-0.21 -0.17 
+0.33 +0.21 
-0.26 -0.15 
+0.52 +0.28 
-0.25 -0.17 
+0.30 +0.37 
-0.28 -0.32 
+0.62 +0.61 
-1.09 -0.58 
+0.45 +0.27 
-0.43 -0.16 



(0.000- 
(0.250- 
(0.500- 
(0.750- 



0.250 
0.500 
0.750 
1.000 



o.68 ±0.17 18:12 1838 
i.io ± o.2i I8i 3 7 18il 

1.52 ± 0.22 18 ?? 18^ 
3.29 ± 0.32 18f 18:1 



i.io ±0.20 181 18:2 

i.63 ±o.25 18:3 18:23 
2.43 ± o.27 1™ 183* 

5.10 ± o.38 l8:i 183 



Table 4: The differential cross sections da/dx® BS for channel (1) as function of x® BS 
for the kinematic region E J T etl > 7 GeV, E 3 T et2 > 6 GeV, as given in Fig. 6, and for the 
kinematic region E J T etl > 6 GeV, E J T et2 > 5 GeV. The 

x obs ran g e j s gi ven i n brackets. The 
quoted cross sections correspond to the centres of the corresponding bins. The first error 
is statistical, the second is systematic and the third one is the energy scale uncertainty. 
Overall normalisation uncertainties due to luminosity measurement (±1.4%) and to D* 
and D° decay branching ratios (±3.7%) are not included in the systematic errors. 
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the total charm dijet cross section is i^g^o. For the x° BS differential cross sections they 
were added in quadrature to the statistical errors and are indicated as the outer error 
bars in Fig. 6. The energy scale uncertainty is shown as the shaded bands. Table 4 lists 
separately the statistical, systematic and energy scale uncertainties. 



7.2 Comparison with Theoretical Predictions 

In Fig. 6(a) the da/dx® BS distributions of the HERWIG MC simulation, normalised to the 
data, are shown for the LO-direct and LO-resolved contributions as well as for their sum. 
The fractions of each contribution was taken from the MC simulation. There is a peak 
in the data at high values of x® BS , consistent with a large contribution from LO-direct 
photon processes. However, there is also a substantial tail to low x® BS values, which 
is not described by the LO-direct MC. Hence a LO-resolved component is required. In 
the LO-resolved MC histogram, the dominant contribution from photon charm excitation 
(lightly hatched) is distinguished from that of other LO-resolved photon processes (densely 
hatched). The contribution of 6-quarks to D* production was taken into account in the 
MC sample as in the inclusive D* analysis. It is about 10% and approximately constant 
with x® BS . The MC distributions, where the LO-resolved and LO-direct contributions are 
allowed to vary independently, were fitted to the data. The data require a LO-resolved 
contribution of 45 ±5 (stat.)%. This value is consistent with the LO HERWIG prediction 
of 37%. The charm excitation contribution to the LO-resolved photon process in the 
HERWIG MC is 93%. 

A comparison of the data with a NLO calculation for a charm dijet sample was performed 
using the massive charm approach [|J. This calculation does not have an explicit charm 
excitation component, since charm is not treated as an active flavour in the photon struc- 
ture function. The x® BS distribution at the parton level was estimated by applying the 
KTCLUS jet finder to the two or three partons produced in this NLO calculation |3(| for 
the kinematic region of our D* and associated dijet analysis. Here e = 0.02 was used and 



m± = \Jm% + (p±), where {p\) is the average p\ of the two charm quarks. The result of 
this calculation (full histogram) is compared to the data in Fig. 6b. To minimise migration 
effects due to hadronisation from high x° BS , the data are given in wider bins compared to 
Fig. 6a. It can be seen that the NLO massive charm calculation Q produces a tail towards 
low x® BS values similar to the light parton jet case fl4i~ |. However, there is a significant ex- 



cess in the data over this NLO prediction. From MC studies we estimate that ~ 6% of the 
highest x° BS bin (0.75 < x° BS < 1.0) can migrate to the lower bins due to hadronisation 
effects. An effect of this size cannot explain the measured low x° BS cross section. Using 
fin = 0.5m^ and m c = 1.2 GeV in the calculation (dashed histogram) yields a higher 
x° BS tail, which is still below the data. With these parameters the cross section near 
x® BS =l is above the data. Applying an intrinsic transverse momentum (h?) = 1 GeV 2 
(section 6.1) increases the predicted cross sections in the two central x° BS bins. However 
the predicted cross sections are still below the measurement. 

The conclusions drawn above are the same when: a) the hadron level jet cuts >6 GeV, 
E?t >5 GeV (Table 4) were used; b) a cone jet algorithm |EJ was applied instead of 
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a cluster algorithm; c) PYTHIA MC was used instead of HERWIG; d) the jet trigger 
described in section 2 was used instead of the nominal one. 

The extent to which a NLO calculation of charm in association with dijets may describe 
the da/dx® BS distribution must await further theoretical developments. In particular 
additional contributions arising from a photon structure function require massless charm 
NLO predictions for da/dx° BS . 



8 Summary and Conclusions 

The integrated and differential inclusive photoproduced D* ± cross sections in ep colli- 
sions at HERA have been measured with the ZEUS detector in the kinematic region 
Q 2 < lGeV 2 , 130 < W < 280 GeV, pf > 2GeV and -1.5 < if* < 1.5. The 
cross section a ep _> D *± x =18.9 ± 1.2 (stot.)lj;g (syst.)nb was measured using the chan- 
nel D* + — > D° -Kg — > (K~ 7r + ) -Kg. A second D* decay channel has been studied, 
D* + — > D° Kg — > (K~ k + k + k~) 7rJ, and good agreement with the Kk: channel has 
been found in the region of overlap (pj > 4 GeV). The results are compared with massive 
and massless charm scheme QCD NLO predictions. The NLO calculations are generally 
below the measured cross sections, in particular in the forward direction. The results are 
sensitive to the parton density parametrisation of the photon used to calculate the cross 
section in the massless charm scheme. 

A sample of inclusive dijet events with an associated D* meson has been used to measure 
the cross section da/dx° BS in the range 130 < W < 280 GeV and Q 2 < 1 GeV 2 . The 
jets were reconstructed with the KTCLUS algorithm, requiring Irf^l < 2.4 and at least 
one D* in the range — 1.5 < i] D * < 1.5 and p®* > 3 GeV. Cross sections are given for the 
kinematical regions E° T etl > 7 GeV, E° T et2 > 6 GeV and E J T etl > 6 GeV, E J T et2 > 5 GeV. A 
peak at high values of x° BS is seen, in agreement with the expectation for direct photon 
processes. A large cross section is also measured at low x® BS , where resolved processes 
are expected to contribute significantly. A comparison of the x® BS distribution to MC 
simulations yields a contribution to the cross section of about 45% from LO-resolved 
photon processes and indicates the existence of charm excitation in the photon parton 
density. The data at x° BS < 0.75 are higher than a NLO massive charm calculation at 
the parton level. 
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Figure 1: AM distributions (a) for the (i-T 77)715 channel with p^* > 2 GeV and (b) for 
the (i^7T7T7r)7r5 channel with p±* > 4 GeV. The full dots are right charge combinations 
from the D° signal region (1.80 — 1.92 GeV). The dashed histograms are wrong charge 
combinations from the D° region for the (K 71)11$ channel and side bands combinations 
(see text) for the (i^7T7r7r)7T5 channel. The full lines are the results of fits to a sum of a 
Gaussian and the functional form A ■ (AM — m n ) B . The insets in (a) and (b) are the 
M(Kir) and M(Kmrrr) distributions from combinations having 143 < AM < 148 MeV. 
The dashed histograms are wrong charge and side bands combinations, respectively. 
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Figure 2: The differential cross section da /dp®* for D* photoproduction, Q 2 < 1 GeV 2 , 
in the kinematic region 130 < W < 280 GeV and -1.5 < r] D * < 1.5 for the (Kir)n s (full 
dots) and (Kmnr)Ti s (open dots) channels. The (Kn)iis points are drawn at the positions 
of the average values of an exponential fit in each bin. The (Ktttt^tts points are offset 
for clarity. The inner part of the error bars shows the statistical error, while the outer one 
shows the statistical and systematic errors added in quadrature. The predictions of NLO 
perturbative QCD calculations are given by the dash-dotted, dashed and dotted curves 
for the massive charm approach J|] and by the full upper (lower) curve for the massless 
charm approach calculation of ref. || (ref. 0), with the parameters described in section 6. 
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Figure 3: Differential cross sections da/drj D * for D* photoproduction, Q 2 < 1 GeV 2 , 
in the kinematic region 130 < W < 280 GeV and a) pf* > 2 GeV; b) > 3 GeV ! 
c) p^* > 4 GeV; d) > 6 GeV. The (Ktc)tts points are drawn at the centres of 
the corresponding bins. The (Khix'k)'k g points are offset for clarity. The inner part of 
the error bars shows the statistical error, while the outer one shows the statistical and 
systematic errors added in quadrature. The curves correspond to the same predictions of 
NLO perturbative QCD calculations as in Fig. 2. 
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Figure 4: Differential cross sections da/dr] D * for .D* photoproduction, Q 2 < 1 GeV 2 , in 
the kinematic region 130 < W < 280 GeV for the (Kn)iTs channel with (a-b)p^* > 3 GeV 
and (c-d)p^* > 4 GeV. The points are drawn at the centres of the corresponding bins. 
The inner part of the error bars shows the statistical error, while the outer one shows 
the statistical and systematic errors added in quadrature. The curves are the predictions 
of the massless charm NLO of ref. || (a,c) and ref. (b,d) with various photon parton 
density parametrisations. 
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Figure 5: Uncorrected transverse energy flow with respect to the jet axis for dijet events 
containing a D* in a kinematic region given in the text and for E^\ al > 4 GeV. The jets 
are defined using the KTCLUS jet algorithm. The distributions are given in three regions 
of rf et separately for direct (x° BS > 0.75) and resolved (rc^ BS < 0.75) photon events. 
The data (dots) are compared to expectations of the HERWIG MC (full histogram) and 
LO-direct only (dotted histogram). The error bars represent the statistical uncertainty 
only. 
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Figure 6: The differential cross section da/dx° for dijets with an associated D* meson 
with pf* > 3 GeV, -1.5 < rf* < 1.5 in the kinematic range 130 < W < 280 GeV, 
Q 2 < 1 GeV 2 , \r] jet \ < 2.4, E J T etl > 7 GeV and E J T et2 > 6 GeV. The KTCLUS algorithm is 
used for the jet definition. The points are drawn at the centres of the corresponding bins. 
The inner part of the error bars shows the statistical uncertainties. The outer part is the 
statistical and systematic errors added in quadrature. The energy scale uncertainty is 
given separately by the shaded bands. In (a) the experimental data (dots) are compared 
to the expectations of the HERWIG simulation, normalised to the data, for LO-direct 
(right hatched), LO-resolved (left hatched), LO-resolved without charm excitation (dense 
hatched) and the sum of LO-direct and LO-resolved photon contribution (full histogram). 
In (b) the data are compared with a parton level NLO massive charm calculation [|J with 
the parameters described in section 7.2. 
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